Theoretical Physics Institute 
University of Minnesota 



TPI-MINN-95/15-T 
UMN-TH-1345-95 
hep-ph/9505289 



Recent Progress in the Heavy Quark Theory 

Talk at the V International Symposium on Particles, Strings and Cosmology 

-PASCOS- 

Johns Hopkins University, Baltimore, March 22 - 25, 1995 



M. Shifman 

Theoretical Physics Institute, Univ. of Minnesota, Minneapolis, MN 55455 

Abstract 



Some recent developments in the heavy quark theory are briefly reviewed. The 
main emphasis is put on interrelation between HQET and OPE. The notion of 
duality and deviations from duality are discussed in detail. 



1 Instead of Introduction 



QCD is a messy theory in the strong coupling regime and with no obvious expansion 
parameter. If we could choose for Nature we would definitely pick up something more 
elegant, preferably in less dimensions, with more symmetry and so on. Alas, the 
choice is not ours. This is the theory of matter, and we have to learn how to deal 
with it as it is. Although there has been no decisive breakthrough in QCD in the 
recent years, a steady progress takes place in different directions. My topic today is 
the progress in the heavy quark theory. 

The reasons why the heavy quark is the best probe of the strong interactions (out 
of a rather scarce supply we have at the moment) are well known. The road which 
led people to the understanding of this fact was not straightforward, though. The 
history of development of the pertinent ideas and formalisms is briefly summarized 
in my recent talk and in excellent reviews ||, and I will save time, sparing the 
audience the n-th repetition of (interesting and important) topics that have become 
routine. Instead, I will concentrate on some very recent results. 

The development of the theory goes in the direction of new applications on the 
one hand, and a deeper understanding of foundations, on the other. Among the for- 
mer let me mention a few which might have a significant impact on phenomenology. 
The first example of this type is the issue of the higher-order perturbative correc- 
tions in the inclusive semileptonic width of B mesons where a noticeable progress 
has been reported |5|. Another example is a new calculation |3j of the power 
corrections to the b —>■ c transition form factor at zero recoil. 

On the theoretical side, we are witnessing how the focus is shifting away from 
building particular 1/rriQ expansions towards more general aspects. Surprising 
though it is, it was not fully realized that HQET is merely a special version 
of Wilson's operator product expansion (OPE) ||. As this fact gains recognition 
more and more elements of Wilson's OPE are being adapted in the HQET environ- 
ment leading to powerful, and sometimes unexpected, consequences. For instance, a 
conspiracy between the infrared renormalons and condensates is common knowledge 
in all problems where the operator product expansion can be used ||. The issue of 



the infrared renormalons in the context of HQET was raised in Refs. [10, 11|. It 



took some effort to reveal peculiarities of HQET and to figure out how concretely 
the same phenomenon of conspiracy works for heavy quark expansions |10, [12] (see 
also |I3| ). The research did not stop at this point, however. 

As was already mentioned, practical needs of phenomenology require inclusion of 
higher-order perturbative corrections. In such problems as the weak inclusive decays 
of the heavy flavors complete calculation of even 0{a 2 s ) terms is a notoriously difficult 
task. People tried playing with different approximate approaches; the so called l b 
graph dominance' prescription due to Brodsky, Lepage and Mackenzie (BLM) JTJj] 
is one of the most popular. Originally the BLM hypothesis was engineered as a 
scale setting procedure in the 0(a s ) terms allowing one to minimize the coefficients 
of the 0{a 2 s ) terms without their complete calculation. This idea was married with 
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the heavy quark theory M. At the next stage the renormalon studies and the 
BLM calculations fuse giving rise to an approximation which is sometimes called 
'resummation of the renormalon chains' [T5|. Using an extended BLM hypothesis 



one isolates and sums up an infinite sequence of graphs which possibly gives an 
estimate of higher-order effects. The results obtained so far || |T3J are only the 
beginning of the story. The extended BLM resummation procedure has not yet 
been properly formulated in the context of OPE. Terms of a very high order in 
the 'renormalon chain' come from large distances and are not calculable in this 
way. The large-distance contribution should be isolated and treated separately, 
as a condensate. The separation is readily done in Euclidean calculations but is 
technically hard to implement in the problems with the Minkowski kinematics, like 
in the inclusive B decays. Until this is done the resummation of the renormalon 
chains is hardly more than an academic exercise and of a limited practical value. 
Still, this is obviously a new promising direction deserving further efforts that will, 
perhaps, eventually find its place in the OPE-based theory. 

First attempts of the instanton calculations in the inclusive heavy flavor decays 
have been reported recently fll6] . The instanton contribution manifests itself in 
the coefficient functions, as a deviation from the so-called practical version of OPE 
WJ\ . This is the only known example of this type which, unfortunately, is very 
poorly controllable numerically. At the present stage one can speak, rather, about 
theoretical aspects of the instanton contribution as a laboratory or a testing ground 
for different effects not seen within the practical version of OPE. The topic is not 
mature enough to warrant a detailed discussion at this Symposium - the existing 
publications leave more questions than answers and should be, rather, viewed as an 
invitation for further work. I will comment, though, on some general aspects of the 
instanton calculations later on. 

The ideas people play with now existed in QCD, in this or that form, for years. 
The heavy quark theory, being one of a few branches of QCD still in active growth, 
gave a new life to them. Being unable to dwell on all these developments I have 
picked up two samples from the current flow: one problem is applied, the other is 
more fundamental. 



2 Deviations from the Heavy Quark Symmetry at 
Zero Recoil and Determination of V c b 

The discussion below will hopefully show how fruitful the intrusions of the OPE- 
based methods are (familiar from the QCD sum rules) for heavy quark theory. 

Needless to say determining V c b is one of the most important current tasks in 
particle physics. Since isolated quarks do not exist it is impossible to measure V c b 
directly, without heavily exploiting a theoretical kitchen. Strictly speaking, there are 
two different kitchens; one of them deals with extracting V c b from experimental data 
through analysis of the exclusive B — > D*lv decays in the so called small velocity 
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limit (slow D*'s). Extrapolation of the amplitude to the point of zero recoil yields 
| Vcb | Fb^d* (zero recoil), where Fb^d* is an effective B — > D* transition form factor. 
In the small velocity (SV) limit this form factor is close to unity as a consequence of 



the heavy quark symmetry [18, [T^, 2C[; deviations from unity are quadratic in the 



inverse heavy quark mass, 1/mf c [JOI ETJ. Our task is to calculate these deviations. 



To predict (Fb-,d* — 1) at zero recoil we derive a sum rule for the transitions 
B — > D* and B —>■ vector excitations generated by the axial-vector current, = 
&7 M 7 5 c. If the momentum carried by the lepton pair is denoted by q, the zero recoil 
point is achieved if q = 0. To obtain the sum rule we consider the T product 

V = if rf 4 xe-^^-(i?|T{4(x)^(0)}|i?) (1) 

where the hadronic tensor h^ v can be systematically expanded in Aqcd/^6,c- For 
our purposes it is sufficient to keep the terms quadratic in this parameter and to 
consider only one out of five possible kinematical structures, namely hi, the only 
structure surviving for the spatial components of the axial-vector current, see e.g. 



22, 23 



Next we use the standard technology of the QCD sum rule approach. Let us 
define 

e = M B - M D * — q — AM - q . (2) 

If e is positive we sit right on the cut. The imaginary part of the amplitude ([!]) 
is the sum of the form factors squared (taken at zero recoil). The sum runs over 
all possible intermediate states, D* and excitations. We want to know the first 
term in the sum, \Fb-,d*\ 2 - Alas, the present-day QCD does not allow us to make 
calculations directly in this domain. 

On the other hand, if e is negative we are below the cut, in the Euclidean domain. 
Here the amplitude ([[]) can be calculated as an expansion in l/m& |C provided that 
|e| Aqcd- To get a well-defined expansion in l/m& )C we must simultaneously 
assume that e <C m h c . 

The non-perturbative corrections we are interested in are due to the fact that 
both, the c quark propagator connecting the points and x in Eq. (JT|) and the 
external b quark lines, are not in the empty space but are, rather, submerged into 
a soft-gluon medium, a light cloud of the B meson. Two parameters characterizing 
the properties of this soft medium are relevant for our analysis. A chromomagnetic 
parameter 

measures the correlation between the spin of the b quark inside B and the chro- 
momagnetic field B created by the light cloud. The second parameter is = 
(2MB)~ 1 (B\b(iD) 2 b\B) measuring the average spatial momentum squared of the 
b quark. Both parameters are proportional to Aq CD . That's all we need for the 
leading non-perturbative term. 
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If the amplitude (|I|) is considered in the Euclidean domain far below the cut (i.e. 
— e ^> Aq CD ) the distance between the points and x is short and we can expand 
hi in AQ CD /m^ c . Actually, the whole amplitude contains more information than we 
need; the sum rule sought for is obtained by considering the coefficient in front of 
1/e in hi. In this way we arrive at the following prediction: 




(4) 

where the sum on the left-hand side runs over excited states with the appropriate 
quantum numbers, up to excitation energies ~ e. (In other words, e plays the role 
of the normalization point. Higher excited states are dual to the graphs with the 
perturbative hard gluon in the intermediate state are neglected together with the 
latter). All form factors in Eq. (|]) are taken at the point of zero recoil. 

Let us now transfer the contribution of the excited states to the right-hand side 
and account for the fact p4|, || that > jjlq. Then we get a lower bound on the 



deviation of F B _+ D * from unity, 



V,-F B ^ D .>^ (5) 



Here we included the perturbative one-loop correction [5R [L9| so that 1 — > rjA 



a s fmh + m c , mu 8\ 
Va = 1 + — — -log—-- ^ 0.975. 6 

7r Krrib — m c m c 3/ 

Using the known value of \Xq and m c =1.3 GeV we conclude that Fb->d* < 0.94. 

Including the — /j,q term and the contribution from the excited states lowers 
the prediction for Fb^d* making deviation from unity more pronounced. If ^ is 



taken from the QCD sum rule calculation p6| the estimate of F B ^ D * is reduced 
to 0.92. As far as the excited states are concerned a rough estimate of the Dn 
intermediate state can be given H implying that 



0.89 ±0.03. (7) 



The error bars here reflect only the uncertainty in the excited states. The parameters 
fj^, fj,Q, m c and t)a have their own error bars which I can not discuss here due to 
time/space limitations. Some relevant numerical estimates can be found in Ref. 



The corrections 0(l/ml c ) to the form factors at zero recoil have been discussed 
previously |28|, |29| within an HQET expansion. In this version, instead of the excited 



state contribution, one deals with certain non-local correlation functions which are 
basically unknown. The advantage of the sum rule approach presented above seems 
obvious: here we deal with only local operators; all non-locality is hidden in the 
excited state contribution which (i) has definite sign; (ii) its magnitude is known to 
be relatively small numerically; (iii) it can be relatively reliably estimated in models. 
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3 OPE and Deviations from Duality 



Many applications of the heavy quark expansions are based, apart from the ex- 
pansion per se, on duality. All problems associated with the inclusive heavy flavor 
decays belong to this class. Needless to say they are important. As usual in QCD, 
a legitimate formulation of the expansion can be given in the Euclidean domain. 
The physically measurable quantities (e.g. the lepton spectrum in the inclusive 
B — > lvX u decays) belong to the Minkowski domain. We want to extract conse- 
quences of the Euclidean expansion in the Minkowski domain. A bridge between 
these two is provided by duality, a concept which is very frequently used but is very 
poorly studied. The notion is not even clearly defined in the literature, so that 
different people ascribe to it different meaning. 

Now I will dwell on some general features of the operator product expansion 
in QCD with the intent to at least initiate discussion of this God-forgotten topic, 
deviations from duality. Many of the arguments presented below are part of ongoing 



research 130 



To establish an appropriate setting and introduce necessary terminology let me 
start from a classic example, the total cross section of e + e~ annihilation to hadrons. 
More specifically, assume, for simplicity, that only u and d quarks exist, they are 
massless, and the "electromagnetic current" has the form 

U= \{^l^-h^d). (8) 

Inclusion of other quarks or the quark masses is a technical issue which does not 
affect the essence of the construction. Also, the particular choice of the current is 
not of importance. 

The central object of our analysis is the T product of two currents 

= i J e^d 4 xT{j,(x)j + (0)} = (q,q u - q 2 g, u ) T. (9) 

For large Euclidean q 2 (q 2 < 0) 

T = Y / C n (q; f i)OM (10) 

where the normalization point fi is indicated explicitly. (Strictly speaking, one 
should deal with q 2 derivatives of Eq. (fL0|)). In what follows we will also use the 
notation Q 2 = —q 2 . The sum in Eq. ([TDD runs over all possible Lorentz and gauge 
invariant local operators built from the gluon and quark fields. The operator of 
lowest (zero) dimension is the unit operator I, then comes the gluon condensate G 2 , 
of dimension four; the four-quark condensate represents the example of dimension- 
six operators (for further details see Ref. PT]). 

The coefficient functions C n (absorbing the short- distance contribution) contain 
both, perturbative and non-perturbative contributions of different types. It is con- 
venient to start from the unit operator whose coefficient accounts for the entire 
perturbation theory. 
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Consider a typical Feynman graph for the vacuum expectation value of de- 
picted on Fig. 1. Assume that the momentum flowing through the graph is Eu- 
clidean and large, Q — > oo. All subgraphs in this graph are assumed to be renor- 
malized at Q. The final result, being expressed in terms of the running coupling 
constant a s (Q) is finite. In the standard calculation of the diagram of Fig. 1 all 
virtual momenta saturating the loop integrals scale with Q as the first power of 
Q. For large Q we use the perturbative expressions for the quark and gluon Green 
functions and one might say, naively, that we arrive at the standard perturbation 
theory for the coefficient of the unit operator, 

C7 = (11) 
i 

where ai are numerical coefficients. As a matter of fact, the above expression is 
not correct theoretically. One should not forget the following: in doing the loop 
integrations in Cj we must discard the domain of virtual momenta below //, by 
definition of C(p). Subtracting this domain from the perturbative loop integrals we 
introduce in Cj power corrections of the type (ji/Q) n by hand. (If it is possible to 
choose p sufficiently small these corrections may be insignificant numerically and 
can be omitted; such a situation will be realized in what is called a practical version 
of OPE, see below). Needless to say that the parameter \i is auxiliary - the final 
result for (T) must be p independent. 

At any finite order the perturbative contribution, eq. (|TTD , is well-defined. At 
the same time, if the number of loops becomes of order l/a s and each loop carries 
roughly speaking one and the same momentum the total momentum Q is shared 
between many lines. This might be called a semi-hard contribution. One can vi- 



sualize it in terms of, say, 'direct instantons' f32|| . Integration over the sizes p of 
the direct instantons is implied; the integral is saturated at p ~ Q~ l . It should be 
emphasized that the correspondence between the Feynman graphs with l/a s lines 
and the semi-hard contributions described above is qualitative. Still, if we rely on 
this correspondence we may say that the characteristic loop momentum scales like 
Q/ InQ. It is quite clear that the contribution of this type has the form 

AC7 ~ exp (-C/a s (Q)) ~ (12) 

where C is some positive constant and the exponent 7 need not be integer. Typically, 
the numerical value of 7 is large. In this way we get a non-perturbative contribution 
to the coefficient function Cj, a non-condensate non-perturbative term. Such terms 
may or may not be important numerically depending on a particular value of Q 
under consideration. (They are disregarded in the practical version of OPE which, 
thus, applies to the values of Q above a critical point). 

This is not the only source of the non-perturbative terms, however. Consider a 
situation in which one, or a few lines in the diagram of Fig. 1 are in a special regime 
- their momenta are soft in the sense that the characteristic off-shellness does not 
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scale with Q. These lines are marked by crosses on Fig. 2. Graphically one can cut 
these lines and treat the product as a local operator. In this way we arrive to the 
sub-leading operators O n , n ^ 0, in the expansion flTOp. For instance, cutting the 
gluon line we get G 2 . The remaining part of the graph where the loop momenta 
scale as Q will yield the coefficient function in front of the corresponding operator 
O n in the form of a series in a s (Q). Taking the matrix element of the sub-leading 
operator produces the condensate non-perturbative correction. (Analogously to Cj, 
the coefficients C n also receive, generally speaking, a semi-hard non-perturbative 
contribution of the type presented in eq. ([121) . By the same token, in calculating 
the perturbative part of C n one must discard the domain of virtual momenta k < p.). 

Finally, it is conceivable that there are no hard or semi-hard lines at all. The 
external momentum Q is transferred from the initial to the final vertex through 
a large number of quanta (growing like a power of Q) and none of the quanta 
carries momentum scaling with Q (Fig. 3). Of course, in this situation one can 
not speak of individual quanta, one should rather use the language of typical field 
fluctuations transferring the momentum Q without having any Fourier components 
with frequencies of order Q. This type of contribution is not seen in OPE truncated 
at any finite order. It is reflected in the rate of divergence of high orders in the 
power series. It is intuitively clear that the corresponding part of the correlation 
function must be exponentially small, ~ exp(— Q). A transparent example is again 
provided by instantons. This time one has to fix the size of the instanton p by 
hand, p ~ po- Then this contribution is C(exp(— Qpo))- The relation between 
exp(— Qpo) piece and the high-order terms of the power (condensate) series is akin 
the connection between exp(— l/a s ) terms and I ~ orders in the perturbative 
expansion. Indeed, proceeding to higher- dimension operators in OPE we force more 
and more lines to become soft, and the limit of only soft lines is correlated with the 
operators built from a large number of the gluon and quark fields. 

I hasten to make the following reservation. It may happen (and actually hap- 
pens) that some of the loop integrations are not saturated at Q or p but are rather 
logarithmic. Logarithms ln(Q/p) occurring in this way are associated either with the 
running coupling constant a s (Q) or with the anomalous dimensions of the operators 
at hand. They will play an important role in the analysis following below. 

Summarizing, Wilson's OPE ([H]) yields (T) in the deep Euclidean domain as 
an expansion in different parameters. Purely logarithmic terms (h\.Q)~~ l are due 
to ordinary perturbation theory. Terms of the type (Q 2 )~ k (\nQ)~' y reflect higher- 
dimension operators and semi-hard contributions ('direct instantons'). In the former 
case the values of k are integer, the latter case may produce non-integer values of 
k. In the practical version of OPE we neglect the semi-hard contributions calculat- 
ing the coefficient functions perturbatively. All non-perturbative terms come from 
condensates within this approximation. 

The perturbative calculation of the coefficient function yields a series in a s . Ev- 
erybody knows that the series is factorially divergent. First, the number of graphs 
grows factorially with the number of loops. This factorial behavior can be corre- 
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lated with the instanton-antiintstanton contribution p3| . Second, there exist specific 



graphs which are factorially large by themselves, the renormalons [34]. As I have 



already mentioned the latter received much attention lately, in connection with var- 
ious applications of HQET. In the framework of the operator product expansion the 
infrared renormalons per se play no role since they are completely hidden in the con- 
densate corrections. By construction, in calculating the coefficient functions we cut 
off the infrared contribution at fi, never come close to the Landau pole and, hence, 
the coefficient functions, calculated properly, are free from the infrared renormalons. 
(The ultraviolet renormalons are not absorbed in the condensate corrections, but 
they are harmless anyway.) 

Passing to the next level, from the perturbative series to the condensate power 
series, we face the same problem since this power series is also asymptotic and is 
truncated in practical calculations (only those operators whose dimension is smaller 
than some number are retained). The uncertainty introduced in this way in the 
Euclidean calculation is expected to be less than the last term retained. As has 
been already mentioned the effect of the high order 'tail' of the power terms is 
exponential. The exponential terms are conceptually connected with the divergence 
of OPE for high-dimension operators and do not show up in the truncated OPE. 
One should proceed, however, with extreme caution estimating the uncertainty due 
to the high order power terms in the Minkowski domain (see below). The existing 
theory tells us nothing about when the exponential terms become negligibly small. 
At this point we have to rely on indirect methods and phenomenological information 
existing in this or that channel. 

Let us discuss a serious theoretical question - whether it is legitimate at all to 
retain in the expansion of the correlation functions logarithmic, power and expo- 
nential terms simultaneously. Indeed, formally any l/(lnQ) 1 term is parametrically 
larger that, say, l/Q 4 , and, moreover, any power term is parametrically larger than 
exponential ones. If so, the theoretical uncertainty due to truncation of the logarith- 
mic series - and we, of course, have to truncate it at some finite order - is formally 
larger than even the lowest-order power correction. 

There are two arguments which seemingly justify the procedure of retaining 
power terms simultaneously with the truncated perturbative series. First, in many 
instances power terms have distinct origin, and the corresponding effects do not show 
up in perturbation theory. This is valid, for instance, in all cases where the chiral 
quark condensate is involved. Another example is due to four-fermion operators 
responsible for the Pauli interference in the non-leptonic weak decays of heavy flavors 
p5| . Second, even if the power terms are not different in their structure from those 
occurring in perturbation theory, there is a strong numeric enhancement of the 
power terms inherent to QCD whose origin is not completely clear at the moment. 
Thus, estimating the uncertainty due to infrared renormalons we get Q~ A with the 
coefficient much smaller than that appearing in Wilsonian OPE. 
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3.1 Practical version of OPE 



In analytic QCD it is rather difficult to carry out in practice the consistent Wilsonian 
separation procedure outlined above. Therefore, in applications, one usually settles 
for what is called the 'practical version of OPE'. In the vast majority of applications 
the practical version works very well numerically ||31|| : notable exceptions are also 
known, however IKIJ . 



What is the essence of the practical version of OPE? Let us elucidate this question 
using the example of the previous section. Assume that we want to calculate the 
vacuum expectation value of T (to be denoted by II) limiting ourselves to some finite 
number of loops. For pedagogical purposes let us consider the two-loop graphs of 
Fig. 1. The corresponding expression for II has a generic form 

II(Q 2 ) = <T> = a s J dk 2 ^F(k 2 /Q 2 ) (13) 

where k is the virtual gluon momentum (the Wick rotation to the Euclidean space 
is implied) and F(k 2 /Q 2 ) is a dimensionless function describing the fermion lines in 
the graphs of Fig. 1 and emerging after integration over the angle variables of the 
vector k. Of course, at k 2 < the genuine expression for II has nothing to do with 
the integral (|T3| ) simply because the gluon Green function has nothing to do with 
1/k 2 . Therefore, in the genuine OPE one should cut off the integral from below at 
/ip excluding the domain k 2 < /i 2 , from the perturbative calculation. This domain 
must be referred to the vacuum matrix elements of different gluon condensates which 
certainly can not be calculated perturbatively. These gluon operators are normalized 
at yUp, by construction, and their coefficients are given, to the leading order in a s , 
by the k 2 derivatives of the function F(k 2 /Q 2 ) at k 2 = 0. The Taylor expansion of 
this function starts from k 2 |3l|], generating the coefficient of G 2 ^ the next term of 
the Taylor expansion of F(k 2 /Q 2 ) will give rise to the coefficient of (D Q G Mt/ ) 2 , and 
so on. Then the coefficient function Cj is given by the integral ([13|) with k 2 > /i 2 ,. 

In the practical version of OPE, instead, one does the following. Let us add 
the integral (|13|) in the domain < k 2 < ^ to the coefficient Cj. One should not 
think that this added integral correctly represents the low k 2 part of the graphs 
at hand. The idea is just to add this integral, by definition, and to subtract it. 
Then the coefficient Cj, to order 0(a s ), is given by full perturbative graphs of Fig. 
1. To avoid double counting, however, we subtract the same contribution from the 
condensates. More exactly, we subtract almost the same contribution. Since the 
sum over condensates is truncated at some finite order n , the best we can do is to 
subtract from the condensate part the integral 



a. 



ru 1 n o 1 

/°4E>(W). (14) 

J K n =2 nl 



In this way we subtract from each condensate its 'one-loop perturbative value'. 

It is quite clear that this strategy of converting the genuine Wilson OPE into 
a practical version (i) leads to a loss of factorization of short and large distance 
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contributions inherent to Wilson's OPE; (ii) can not be systematically generalized 
to any number of loops since there is no way to unambiguously define the integrand 
to all orders in the small k 2 domain (iii) only approximately avoids double counting 
due to the necessity of truncating the condensate series at a finite order (the lower 
the order we truncate the larger error is introduced; likewise, the larger value of /ig 
is chosen the larger error is introduced). 

Therefore, the practical version described above is useful (i.e. sufficiently accu- 
rate numerically) only provided /z 2 , can be chosen small enough to ensure that the 
'one-loop perturbative' contributions to the condensates are much smaller than their 
genuine values and, at the same time, a s (/ig)/7r is small enough for the expansion 
to make sense. The existence of such a window is not granted apriori and is a very 
fortunate feature of QCD. The practical version obviously ignores the semihard con- 
tributions to the coefficient functions discussed above. Working within this version 
we automatically limit ourselves to a finite order in a s forbidding any questions 
concerning any aspects of the high order behavior. 



3.2 Consequences of OPE in the Minkowski domain 

Now, let us recall that in the problems at hand we are not interested in the correlation 
functions at the Euclidean values of q 2 per se. The physically observable quantity is 
the imaginary part. The question is what can be predicted about the imaginary part 
from the Euclidean expansion of the type discussed above. In other words, we would 
like to reinterpret the expansion in terms of the imaginary parts it corresponds to. 

If in the Euclidean expansion the notion of the soft and hard lines means small 
and large off-shellness, respectively, one can not extend this classification to the 
Minkowski domain. This is especially clear for the perturbative contribution which, 
by means of the Cutkosky rule, is obtained by putting all quarks and gluons in the 
intermediate state on shell. 

Of course, quarks and gluons are confined and have no mass shell, literally speak- 
ing. The genuine imaginary part is saturated by hadrons, not quarks and gluons. At 
large s (where s = q 2 ) the numerical value of the genuine imaginary part coincides 
with that obtained by calculating with quarks and gluons only approximately. How 
large is the actual deviation between the genuine imaginary part and that inferred 
from the quark-gluon calculation in the Euclidean domain? This is the key question 
the QCD practitioners have to address sooner or later. I will suggest some tentative 
answers. 

Let us first explain what is meant when one speaks about the quark-gluon imag- 
inary part. The starting point is the Euclidean calculation. Consider first the 
perturbative series for Ti(Q 2 ) truncated at some finite order below the critical order 
where the perturbation theory becomes senseless. The successive terms in this series 
have the form 

1 1 , lnlnQ 2 

and — — — — - and so on. 



InQ 2 ' (InQ 2 ) 2 ''"' (InQ 2 )' 
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Now, in each separate term of this series we substitute Q 2 — > — s and then take the 
imaginary part. 

We repeat the procedure with the condensate terms, again truncated at some 
finite order below the point where the asymptotic divergence shows up. In the 
Euclidean domain a generic contribution due to higher dimensional operators is 

2n 7 



(Q 



2\k 



x [(lnQT + - 



where k is an integer number related to the (canonical) dimension of the operator 
at hand. It translates into 

in the imaginary part. Notice that the imaginary parts concentrated at s = (i.e. 
5(s) and its derivatives) do not manifest themselves at large s in the truncated series. 
Notice also that in the absence of the logarithmic factors (i.e. for 7 = 0) we get no 
imaginary part at large s. This feature is specific to the example at hand and does 
not extend, generally speaking, to other problems. 

Summarizing, we do a straightforward analytical continuation, term by term. 
The sum of the imaginary parts obtained in this way will be referred to as the quark- 
gluon cross section. This quantity will serve as a reference quantity in formulating 
the duality relations. When one says that the hadron cross section is dual to the 
quark-gluon one the latter must be calculated by virtue of the procedure described 
above. 

Defining the analytic continuation to the Minkowski domain in this way it is not 
difficult to see that those lines which were far off shell in the Euclidean calculation 
remain hard in the sense that now they are either still far off shell or on shell but 
carry large components of the four-momenta, scaling like ^/s. 

The contributions left aside in the above procedure are related, at least at a 
conceptual level, to high-order tails of the both series, the logarithmic a s series 
and the power (condensate) one. These tails generate negligibly small terms in the 
Euclidean domain, falling off at large Q 2 as a high (and, generically, non-integer) 



•Oj 



power of 1/Q or exponentially, ~ exp(— Qp, 

Both can be visualized through instantons. In the first case we deal with the 
direct instanton contribution integrated over p (the p integration is saturated at 
p ~ Q^ 1 )- In the second case p is to be considered as a fixed parameter not scaling 
with Q~ x . As we will see below, the fixed-size instanton is actually relevant in 
both cases, so we start our discussion assuming that p is fixed. The corresponding 
contribution to II(Q 2 ) was found long ago [[36], |37|. In the Euclidean domain it is 



representable as an integral over an auxiliary parameter over a McDonald function. 
We do not intend here to rely on specific details of the instanton expression; rather, 
we would like to abstract general features inherent to mechanisms of this type. To 
this end we will follow the suggestion of Ref. p8| |. The following observation was 



made in this work: if the condensate terms correspond, in the coordinate space, to 
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singularities of U(x) at x = 0, the instanton-type contributions are connected with 
the singularities of U(x) located at a finite distance from the origin; the distance 
from the origin plays the same role as the instanton radius. For our purposes it is 
reasonable to assume that these singularities have the simplest possible structure, 
namely, 

— -, In (x 2 + p 2 ), (x 2 + p 2 ) In (x 2 + p 2 ), and so on. 

x 2 + p 1 

The Fourier transforms of these expressions have the generic form 

(Qp)- n K n (Qp), n = l,2,... 

where K n is the McDonald function. 

As was expected, at large Euclidean Q 2 the corresponding contribution dies off 
exponentially, 

An = (Qp)~ n K n (Qp) oc {Qp)- n - l ' 2 e- Q P , (15) 

in full accord with intuition regarding transmitting a large momentum through a 
soft field fluctuation. We are interested, however, in the Minkowski values of the 
momentum. The question is what is the impact of the contribution above on Imll. 

Analytically continuing All given in Eq. (|1|) to the Minkowski space (Q — > iQ) 
we arrive at 

ImAn = (-l)-+i|(Vip)-"J n (yip) oc (-^^^(^-""^cos^p-^) , (16) 

where J n is the Bessel function. We see that the imaginary part falls off rather 
slowly, as a modest power of l/y/s, and oscillates. 

This result may or may not seem surprising. It tells us that with no further 
dynamical input, on general grounds alone, one can not rule out large violations 
of duality. If this were the end of the story the operator product expansion would 
be, to a large extent, useless in all problems where one needs to predict quantities 
referring to the Minkowski kinematics, such as decay probabilities, spectra, etc. 

Fortunately, the situation is not as bad as it might seem at first sight. Indeed, 
so far we were discussing transmitting a large momentum through a fixed-size fluc- 
tuation. In QCD the size of the fluctuation is not fixed; rather, we integrate it over 
with a weight function which depends on p. 

The result of this integration depends on the choice of the weight function w(p). 
Leaving aside inessential details we are basically left with two distinct possibilities 
- one will represent the tail of the a s series, the other the tail of the condensate 
series. The first option is related to the small-size direct instantons p ~ l/v^j it 
generates a duality violating contribution falling off as a high power of l/y^. The 
second type of the duality- violating contribution, falling off exponentially with y/s, 
is due to instantons whose size is smeared around some po (i.e. parametrically small 
p are excluded). 

Let us start from the small-size instantons. If yfs is parametrically large and 
p ~ 1/ t/s is parametrically small one can use the dilute instanton gas approximation 
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Then the instanton density d(p) (related to the weight function mentioned 
above) has the form 

d(p) = (Aqcdp) 6 (17) 

where b is the first coefficient in the Gell-Mann-Low function. For three (almost) 
massless quarks 6 = 9. For small p the weight function w(p) has the same form as 

w(p) = (A QCD p) 7 /(p), (18) 

but the exponent 7 depends on dynamical details. Say, if we treat three light 
quarks as massless (i.e. neglect their mechanical masses) then w(p) is proportional 



to the cube of the quark condensate E0[ and 7 = 18. If the strange quark mass is 



considered as a parameter < Aqcd but the mechanical masses of u and d are still 
neglected w(p) is proportional to the square of the quark condensate and m s , and 
and 7 = 16. If all three light quarks are treated as massive 7 = 12. Even following 
the latter (quite unrealistic) treatment we can safely say that 7 is large numerically. 
Let us remind that generically 7 need not be integer. For instance, if one considers 
very small values of p where the c quark can be treated as light (i.e. pm c <C 1) one 
should use 7 = 55/3 in the chiral limit with respect to u, d, s quarks. This example 
is quite academic, though, since at such small values of p the weight function is so 
strongly suppressed that the corresponding contribution can hardly be observable 
in hadronic processes in the foreseeable future. 

Since our purpose is mostly illustrative we will not restrict ourselves to a specific 
value of 7; rather it will be treated as a generic numerical parameter lying somewhere 
in between 10 and 20. 

The function f(p) in Eq. (|18D is a cut off function suppressing large values of 
p. It is equal to unity in the dilute instanton gas approximation. In fact, one does 
not need any cut off if calculations are done at large Euclidean momenta Q 2 . The 
integral is automatically convergent at p ~ 1/Q (see Eq. (|T^)). Since we deal now 
with the imaginary parts that fall off at large yfs only slowly (see Eq. ([TB])), to be 
on the safe side, we introduce a cut off function. The final result will be essentially 
independent of the choice of f(p) provided the cut off is soft enough. For instance, 
f(p) = exp(— p/po) is quite suitable. 

With this information in hands we proceed to smearing ImAII given in Eq. (|16|), 



(ImAn> s oc / (Vip)--J n (VipMp)- (19) 
J p 



where the angle brackets denote smearing and the subscript s indicates the we 
average over the small size fluctuations. Doing the integral we arrive at 

(ImAH). oc (sin =2) 2-r(2)r(2^) x (20) 

plus subleading in terms. This expression assumes that y/sp 3> 7, and 

characteristic values of p saturating the integral (|19|) are of order p ~ l/yfs- The 
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very same expression for the imaginary part could be obtained by integrating the 
Euclidean All over p with the subsequent analytic continuation to the Minkowski 
domain and separation of the imaginary part. 

What lessons can we learn from this simple exercise that describes, at a qual- 
itative level, transmitting large momenta through a small-size fluctuation? First, 
this particular duality- violating contribution to the imaginary part falls of as a high- 
power of 1/a/s at large s. In order to be able to single out the small-size fluctuations 
one must choose -Js > 7/po- Even if p$ 1 is set at 0.5 GeV (in fact, this parameter is 
probably larger |4l[| ) the small-size direct instantons can show up at energies y/s > 
5 GeV. At these energies 



<ImAn> s ~ < (l/2)f (21) 




and this effect is unimportant for all practical purposes. Parameter Aqcd in Eq. 
(pl|) can be hardly pinned down better than up to a factor of two at present but 
this is quite irrelevant - the nature of the s dependence in Eq. ( plj) is such that 
we can safely approximate it by the step function. If ^/s is at least factor of two 
higher than 7Aqcd the result is practically zero; shifting slightly towards lower 
values of a/s we formally observe an almost immediate explosion in Eq. (|2T|). This 
explosion happens, however, outside the domain where one can speak about the 
small-size fluctuations. The conclusion is as follows: immediately above a critical 
point (ImAII) s practically vanishes. Below the critical point no reliable predictions 
exist in the Minkowski domain. 

Let us discuss now another contribution which is formally much smaller than 
that of the small-size instantons, namely the contribution coming from the vacuum 
medium where the size of the fluctuations can not be arbitrarily small. This is 
presumably the dominant component of the vacuum fields responsible for the most 
salient features of the hadronic processes. Being formally exponentially suppressed 
this mechanism at the moment seems to be the leading source of the duality viola- 
tions in all processes of practical interest. To have something particular in mind one 



can think of the instanton liquid model [41], although this model is certainly not 



singled out in the context under discussion. The general features will be inherent 
to any mechanism of this type. 

As was already demonstrated, the contribution of the fixed-size instanton falls 
off very slowly in the Minkowski domain. Certainly, it is not realistic to saturate by 
fixed-size fluctuations; the ensemble of fluctuations in the genuine QCD vacuum is 
characterized by a distribution in the sizes centered near some typical size po- Let 
us assume that this distribution falls off steeply at smaller and larger values of p, 
with a width A. A typical distribution of this type can be parametrized as follows: 

w(p)=Afexp{ /3p} (22) 

P 
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where Af is a normalization constant, 



p is the center of the distribution and A is its width. 

Averaging Im All over this weight function one obtains (for n = 1) 



(Im An); = M2E- 1 J l {V2^[^f3 2 + E 2 - P} 1 / 2 )K 1 {V2^[</p~ 2 + E 2 + B} 1/2 } (24) 

(similar expressions for n > 1 can be derived by differentiating with respect to a; 
the subscript I indicates that we average over the large size fluctuations). If 

E » if (25) 



the smeared imaginary part reduces to 

(Im An), - 2AfE- 1 J l ( s /2^ ^)K 1 (^2Ei ^) . (26) 

It is quite clear that the exponential suppression of (Im An)/ that emerged after 
the smearing is due to the oscillations of the original imaginary part generated by 
the fixed-size fluctuation. If the weight function is narrow (A <C p ) this suppression 
starts at high energies, see Eq. fl25[) . In the limit when A — > with E fixed the 
exponential suppression disappears from Eq. fl2~4]), and we return back to the original 
oscillating imaginary part. Notice also that the exponent at E ^> 4^ is different 
from that one deals with in Euclidean domain [\fE versus Q). Moreover, apart from 
the exponential suppression we observe residual oscillations as well. Although the 
particular results (|25"D, ( f2~6| ) for the duality violating contributions were obtained by 
using instanton-inspired calculations (and are dependent on the choice of the weight 
function) the general features of the results - an exponential suppression starting at 
a point correlated with the width of the smearing function plus residual oscillations 
- are presumably common to all mechanisms and will be confirmed in the future 
solution of QCD (if it is found) and/or phenomenologically. I make a conjecture that 
the exponentially suppressed duality violating contributions have a generic form 

exp(-f(E)), f(E) -> kE u at large E (27) 

where the critical index a < 1. Shortly we will discover the same qualitative picture 
in a totally different context. 

It is worth emphasizing once more that although at academically high energies 
the exponentially suppressed contribution dies off faster than that of the small-size 
instantons in all practical problems the latter is totally negligible once we cross the 
critical point while the former, being a less steeper function, can play a role. 
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3.3 Further explorations of deviations from duality 

Let us consider now a model spectral density suggested in Ref. 0. Although this 
model is not derived as a solution of QCD it satisfies all general properties of QCD 
we are aware of today and seems to be very instructive in studying the duality 
violations. It will be seen that qualitatively it nicely matches with the discussion 
above. 

To begin with we remind the model and its motivations. The model can be 
formulated, in a most straightforward way, for the spectral density associated with 
one heavy and one light quark in the limit when the heavy quark mass tends to 
infinity. Simplifications occur due to the fact that the heavy quark Green function 
then is trivially known. Specifically, let us consider heavy-to-light quark currents 

J s = Qq, J P = Q>y 5 q; 

Ji = \{Js + Jp) = Q\^ + 7b)<Z, -h = \(4 ~ 4>) = + J 5 )Q . (28) 

assuming that m,Q — > oo and m q — > 0. Note that J 2 is not hermitean conjugate to 
J\\ rather J 2 is a chiral partner to j\. Therefore, the correlation function 

n = if e^dxivaclTiJ^x), J 2 (0)}\vac) (29) 

vanishes in perturbation theory if m q is put to zero. 

First of all, we should choose the external momentum k in a most advantageous 
way. If rriQ —>■ 00 it is clear that the optimal reference point lies "slightly" below 
threshold, 

ko = m Q -e, k = 0, 

where e = C(A QC d)- 

Taking the limit niQ —* 00 and neglecting the hard gluon exchanges one can 
rewrite the correlation function 11(e) at positive (Euclidean) values of e in the fol- 
lowing way |JJ: 

j roc rr 

n(e) = -/ e-"dr(q(r)exp{ igA (t)dt}q(0)) (30) 
4 jo J 

where the angle brackets denote the vacuum averaging. The correlation function of 
interest is thus nothing else than the Laplace transform of the Wilson line (in the 
time direction) with the light quark fields at the end points. 
If e ^> Aqcd one can expand Eq. (p0|) in 1/e, 



e z e 4 e° 



(31) 



where Vo is the time component of the (Euclidean) momentum operator. 

From Eq. (|3lD it is perfectly clear that the operator product expansion for 11(e) 
must be asymptotic; otherwise it would define an odd function of e, which is obvious 
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nonsense. Indeed, at positive e we are below the cut, and 11(e) is analytic in e. At 
negative e, however, we sit right on the cut generated by the intermediate physical 
states produced by the currents Ji )2 , and the correlation function II (e) develops an 
imaginary part, a discontinuity across the cut. Qualitatively we have a pretty good 
idea of how Imll looks like. We will return to this issue later on and now observe 
only that if the series in (|31|) is truncated at any finite order the only imaginary 
part we obtain from Eq. ( [31]) is concentrated at e = 0. In accordance with our 
formulation of duality we conclude that at large positive values of E (E = — e), 
E ;§> Aq CD , the physical spectral density is predicted to vanish. 

While the exact spectral density is unknown, of course, a model which should 
be close to it, at least qualitatively, has been suggested in Ref. On general 

grounds we know that at large r the Wilson line (<?(t) exp{/ T igA (t)dt}q(0)) must 
be proportional^ to exp{— At}; moreover, this function should have no singularities 
in the complex plane except on the imaginary axis. At small r it must be expandable 
in powers of r and this expansion must run over only even powers of r (this expansion 
generates Eq. ([H])). The simplest choice satisfying the above criteria is 

(q(r) exp{£ igA o (t)dt}q(0)) = (qq) x — (32) 
This expression for the Wilson line generates a correlation function 11(e), 

n <«> = tk^> 3 par) < 33 > 

where (3 is related to Euler's if) function, 



fc=0 



-1 



x + k 



We pause here to make a digression illustrating the degree to which realistic 
details can be expected to be reproduced by the model above. Accepting Eq. (p2[) 
we made a specific assumption about the behavior of the matrix elements appearing 
in Eq. (|3"T| ) - as dimension of the operators grows their matrix elements grow 
factorially, 

(qCPo) 2n q) - (qq)(A 2 QCD ) n C 2n (2n)\ , (34) 

where C is a numerical constant. Although very little is known for fact for high- 
dimension operators one can confront the dimension-5 operator with information 
that already exists. Indeed, 

(qPoq) = I (qwG%a, u t a q) = \ml (qq) (35) 



1 Warning: in many publications it is erroneously assumed that (q(r) exp{/J" igAo(t)dt}q(0)) oc 
exp(— CV 2 ) which is incompatible with general principles. 
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where Wq is a parameter introduced in JO] 0, and we accounted for the fact that 
Vq is defined in Euclidean. By comparing the first subleading term in the large e 
expansion of Eq. ( |3"3"D with the general expansion ( |5T] ) we establish that (i) the sign 
is correct; (ii) A 2 = mo/8. Phenomenologically m 2 , is close to 0.8 GeV 2 fi2|l , and we 
conclude that in our model 

A » 0.31GeV. 

This number is in the right ballpark; perhaps, a factor of ~ 1.5 smaller than the 
expected phenomenological value of A. 

The correlation function (^) we end up with is a sum of simple poles at 

E = -e = 1,3,5,7,.... 

(in the remainder of this section we put A = 1; thus all energies will be measured 
in these units), all residues are equal in the absolute value and are sign-alternating. 
We remind that \mH(E) is actually the difference between the spectral densities in 
the scalar and pseudoscalar channels and, therefore, need not be positive. 

The above model perfectly matches a picture one expects to get in the multicolor 
QCD, with N c — > oo: two combs of infinitely narrow peaks sitting, back-to back, on 
top of each other. Let us defer for a while the discussion of how the spectral density 
evolves from N c = oo to N c = 3 and study now the implications of this picture for 
duality. 

The spectral density at E 3> A QCD stemming from Eq. (|33|) 

oo 

ImH = -(qq) ^(-l) k 6(E - (1 + 2k)) , (36) 
1 k=0 

taken at its face value, does not vanish. We understand, however, that under the 
circumstances duality should be applied to the average spectral density rather than 
to ImU(E) itself. If we merely average ImU(E) in the interval from E\ to E2 we get 
a sign-alternating result of order 1/(E2 — E{) - i.e. a huge deviation from duality 
which does not even die off as E — > 00 as long as the size of the smearing interval 
is kept fixed! 

If this were the actual situation in QCD any OPE-based expansions in the 
Minkowski domain would be hopeless. For instance, why then bother about cal- 
culating the 1/rriQ corrections to the lepton spectrum in B — > lvX u due to \j? g and 
/i 2 [ 53 1 if uncontrollable deviations may die off slower than l/m 2 ^? Fortunately, at 



large energies QCD with N c = 3 does not look at all like its multicolor limit; the 
limits E — > 00 and N c — > 00 are not commutative. If the number of colors decreases 
from infinity to 3 the spectral density Imll(.E') at large E experiences a dramatic 
evolution. A comb of delta functions is converted into a smooth function. The de- 
tails of the smearing mechanism are foggy but qualitatively one may conjecture the 



2 For a review see the Reprint Volume cited in Ref. plfl ; note that Eq. (2.4) on page 22 of the 
above Volume contains a misprint, — g must be substituted by ig; moreover, a^ v is understood as 
(1/2) b/^- 



following. At N c = oo the correlation functions in QCD are saturated by a master- 
field f44jl , a classical field configuration which does not fluctuate. This freezes the 
scale and makes the calculation in a sense similar to that with the fixed-size instan- 
ton. With N c decreasing, the field configurations saturating the functional integral 
start fluctuating; correspondingly, a distribution of scales emerges. The correlation 
function corresponding to N c = 3 can be obtained, qualitatively, by substituting 
e 



ep in Eq. ([33]) and integrating over p with a weight function similar to (|22|) 
centered at p — 1. In numerical exercises below we will consider the weight function 
(22) as a particular example, with 



a = (3 



A 5 



Af 



2Ki(2/ A 2 ) ' 



where A is the width of the distribution, A = 0(1/N C ). 
Thus, our model spectral density becomes 



/•oo 7]- °° 

ImU(E) = / dpw(p)-(qq) £(-1)** [Ep - (1 + 2k)) 
Jo 2 fc=0 



^fe> £(-i)< 



w 



k=0 



l + 2k s 



(37) 



Analytically it is not difficult to prove that this imaginary part falls off faster 
than any given power of 1/E. Indeed, let us focus on the function 



= £(-1)" 

k=0 



W 



l + 2k 
E 



E 

fc=0 



w 



1 + 4fc 
E 



w 



3 + 4fc 
E 



(38) 



To show that at large E terms 0(1/ E) are absent in cr(E) we represent 



w 



'3 + 4k s 
E 



Then 



-1 oo 
Z k=0 



w 



w 



n + Ak 

1 + 4k\ 



J + w 



w 



5 + 4k 
E 

'5 + 4k s 
E 



0(l/E 2 



OiljE 1 



= w(l/E) + 0(l/E 2 ) 
Moreover, in the next order one uses 



w 



3 + 4k 
E 



w 



l + 4k" 
E 



J + 2w 



5 + 4k 
E 



1 

3 W 



9 + 4k 
E 



J 



+ 0(l/E 3 ). 



an so on. 

Hence, the suppression is exponential, modulated by oscillations. These facts 
are quite general and do not depend on the particular choice of the weight function 
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provided that the weight function falls off faster than any power of p at small p. For 
instance, to prove that an infinite number of oscillations must necessarily take place 
in the imaginary part we write the dispersion relation for 11(e), expand it in 1/e and 
compare the result of this expansion with Eq. (|3TD . The absence of all even powers 
in 1/e requires an infinite number of oscillations in ImIT(i?). 

I was unable to find analytically the type of the exponential suppression. Nu- 
merically it seems that Eq. (|27|) goes through, with the index a depending on the 
width A. (Thus, at A = 1/3 we get a » 0.7). The plots of ImU(E) for A = 0, 
A = 1/3 and A = 1/2 are presented on Fig. 4. We see that the heights of the oscil- 
lations rapidly fall off starting with some boundary energy which, in turn, depends 
on A. Duality is restored with the exponential accuracy at E 1/A. Below this 
boundary value the resonance peaks are conspicuous. Notice also that the length 
of oscillations grows with energy. At A — > the boundary value of energy shifts to 
infinity, and we return back to the infinite comb of the infinitely narrow peaks. 

Theoretically the precise shape of the fall off of \mH(E) at large E is correlated 
with the divergence of the high-order terms in Eq. ([H]). This divergence depends, 
in turn, on the behavior of the matrix elements (qVQ n q) which is not known. In 
our comb-like model (see Eq. ([32])) these matrix elements grow as (2n)\. After 
smearing this formula changes, the result being dependent on the particular choice 
of the smearing weight function. The crucial parameter is nA. As long as this 
parameter is smaller than unity the leading behavior remains intact, (2n)!. For 
nA 3> 1 the regime changes and (2n)! is substituted by A 2n ((2n)!) 2 . We remind 
that A ~ 1/N C . This means that the terms in the operator product expansion ( |3~T|) 
which are subleading in 1/N C are more singular in n, i.e. their divergence in high 
orders of OPE must be stronger. This aspect is not new, though; the very same 



situation takes place in the ordinary perturbation theory [45] 



I am grateful to J. Bagger and A. Falk for inviting me to give this talk at this 
wonderful Symposium. 
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